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Abstract 

Photonic-crystal-based integrated optical systems have been used for a broad range of sensing 
applications with great success. This has been motivated by several advantages such as high 
sensitivity, miniaturization, remote sensing, selectivity and stability. Many photonic crystal sensors 
have been proposed with various fabrication designs that result in improved optical properties. 
In parallel, integrated optical systems are being pursued as a platform for photonic quantum 
information processing using linear optics and Fock states. Here we propose a novel integrated 
Fock state optical sensor architecture that can be used for force, refractive index and possibly local 
temperature detection. In this scheme, two coupled cavities behave as an “effective beam splitter”. 
The sensor works based on fourth order interference (the Hong-Ou-Mandel effect) and requires a 
sequence of single photon pulses and consequently has low pulse power. Changes in the parameter 
to be measured induce variations in the effective beam splitter reflectivity and result in changes to 
the visibility of interference. We demonstrate this generic scheme in coupled L3 photonic crystal 
cavities as an example and find that this system, which only relies on photon coincidence detection 
and does not need any spectral resolution, can estimate forces as small as 10“^ Newtons and can 
measure one part per million change in refractive index using a very low input power of 10“ 

Thus linear optical quantum photonic architectures can achieve comparable sensor performance to 
semiclassical devices. 
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I. INTRODUCTION 


Integrated photonics based on photonic crystal (PhC) structures provides a path to ex¬ 
tremely small optical sensors with applications to biology Bia. chemistry |3] and engi¬ 
neering |3]. PhC devices with various geometries and structures such as hollow core PhC 
hbers [S], ID and 2D waveguides PE] and nano-cavities [81-fT^ have been fabricated and 
used for sensing applications. Among these devices, PhC cavity-based sensors offer impor¬ 
tant advantages over PhC waveguide sensors since they can be made much smaller, thus 
reducing vulnerability from impurities and losses. Moreover, exploiting high Q cavities with 
large mode volume are advantageous for sensors based on refractive index (RI) changes, 
for example in bio-pathogen detection [13], chemical sensing [T3| and single particle detec¬ 
tion [131 . All these schemes make use of second order interference of coherent states of the 
optical held and are thus basically classical phenomenon. 

In parallel to integrated optical sensors, considerable progress has been in using integrated 
optical systems for single photon optical quantum computing using linear optics, so called 
LOQC [IE]. These schemes are enabled by the uniquely quantum-optical phenomenon of 
Hong-Ou-Mandel (HOM) interference [T7]. This opens up a new perspective for optical 
quantum metrology that combines ideas from photonic crystal sensors with linear optical 
quantum information processing using single photon states. We make a hrst step in this 
direction by presenting a scheme that uses HOM interference to make sensors from optical 
pulses prepared in single photon pulses. As HOM interference does not arise for coherent 
optical pulses, our proposal is a true quantum metrology scheme and realizes the gain 
in sensitivity such schemes offer. Furthermore, it opens up the prospect of using LOQC 
protocols to construct more sophisticated quantum metrology protocols that are compatible 
with integrated optical systems. 

Recent demonstrations of cutting edge sensors that exploit quantum mechanics have 
been shown to outperform their classical counterparts in achieving higher sensitivities [lE]^ 
[20] . Many applications, e.g. biological sensing [2T], require low power to preserve delicate 
samples destroyed by: photo-decomposition, photo-thermal effects, and photon pressure for 
example. This requirement is in addition to the usual requirements of high input-output 
gain (responsivity), low noise and high bandwidth. In that regard, weak coherent light offers 
a route to low power sensing. However, the use of weak coherent pulses lowers a sensor’s 
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bandwidth. Consider for example a series of weak coherent pulses with on average one 
photon per pulse, in this case roughly 37% of pulses have no photons at all and 26% have 
more than one photon per pulse. Clearly the ultimate low pulse power limit is achieved 
by single photon pulses with only one photon per pulse. A sensor operating with single 
photon states offers low power suitable for deployment in lab on a chip applications 1221 and 
compatible with attojoule all-optical switching [23] and opto-mechanical devices for strain 
sensors [23| and accelerometers [25] . 

While single photon states are not easy to make there is a very large research effort 
underway driven by their potential application in quantum information processing [16]. For 
our purposes it suffices to note that PhC devices are compatible with a number of quantum 
dot single photon sources [26] and that technological advances in integrated multiplexed 
single photon sources in PhCs are very encouraging [27]. The fundamental quantum nature 
of photons is usually observed through the HOM effect which has now been demonstrated 
in a variety of physical systems such as evanescently coupled optical waveguides [2H] and 
microwave devices [29]. In the HOM effect indistinguishable photons simultaneously arrive 
at each of the two input ports of a 50/50 beam splitter, after which the photons “bunch” 
together so that both photons are either in one output port or the other. Never will you 
observe one photon in both outputs. 

The use of dual Fock states was proposed in 1993 in the context of quantum metrology 
to reduce the uncertainty of phase measurements [30]. In this article, we propose a novel 
scheme for a quantum photonic sensor based on coupled PhC cavities that exploits the 
HOM effect, shown in Fig. The coupled PhC cavities form an “effective beam splitter” 
for two incident photons. The central idea is that a parameter to be estimated, call it "0, 
modulates the coupling between the optical cavities, g. This can be done by changing the 
distance between the cavities through compressing or stretching the dielectric material (e.g. 
for force and strain sensing) or by changing the refractive index of the media between the 
two cavities (e.g. for RI, temperature and single particle sensing). The change in g modihes 
the reflection and transmission of the effective beam splitter which changes the visibility of 
HOM interference. Therefore, by measuring the change in HOM visibility, we can sense the 
variation in g and thus estimate 0. This scheme is independent of transmission/reflection 
spectra normally used for classical cavity-based sensors [31] and neither a dispersive element 
nor spectral resolution for the measurement is required. 
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FIG. 1. Schematic of quantum PhC sensor. Coupled PhC resonators implement an effective beam 
splitter interaction between the input single photon fields resulting in HOM interference effect 
observed in detected output fields. is a measure of the number of coincidences which is 

a function of the time shift between the photons entering the beam splitter. By compressing or 
stretching the distance between optical resonators or through changes in refractive index of the 
medium between the resonators, the coupling between the cavities changes. This results in a change 
in transmission and reflection of the beam splitter and therefore results in a change in the measured 
HOM visibility. 

First, we characterize the proposed sensor in terms of its performance metrics, the respon- 
sivity and minimnm detectable valne for the parameter to be estimated. This characteriza¬ 
tion in terms of the working parameters of the sensor is expressed in a general way, with no 
assnmptions set for the values for the cavity damping rate, cavities coupling strength, PhC 
refractive index, etc. Then, a more specific example is provided by considering this scheme 
with previously reported experimental parameters for GaAs/AlGaAs PhC structures. We 
theoretically predict that such a system can measure one part per million change in refrac¬ 
tive index as well as forces on the order of 10“^ N. These results are not obtained by using 
experimental values specihcally optimized for our scheme. However, the results obtained for 
refractive index and force sensing are promising for integrated on-chip sensing. 


II. HONG-OU-MANDEL SENSOR 

As HOM interference is a uniquely quantum mechanical phenomenon we must necessarily 
proceed with a full quantum description. Consider the double optical resonator scheme 
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composed of two optical cavities with resonance frequency uj, depicted in Fig. The 
optical helds in the two cavity modes are described by the bosonic annihilation operators oi 
and 02 - The interaction picture Hamiltonian is given by 

Hj = hg{a\a2 + aial), (1) 


where g is the effective interaction strength that depends on the parameter -0. We couple this 
system, via the evanescent held, to the input and output channels comprising two optical 
wave-guides. The relation between the respective input and output helds is [32] 




( 2 ) 


where Kj{j = 1 , 2 ) is the damping rate of cavity j. Cavity modes ai{t) and 02 (t) can be 
related to the input modes using the input-output stochastic diherential equations [32l [33] 

. /tj , . , — , . 

—— = -iga 2 {t) - -ai{t) -F V^apinH), 
at 2 

- ^Q 2 (^) + \/^Q 2 ,in(^), (3) 


where the solution to these equations is given in appendix A. To operate this device as 
a sensor we then load the two input ports with single photons and perform coincidence 
detection at the outputs. The upper input port of the beam splitter is loaded with a 
single photon in the state | 1 ^) = / (is^(s)a| j^(s)| 0 ) having pulse shape ^{t) = with 

the normalization condition j dt\^{t)\‘^ = 1. The lower input port is loaded with a single 
photon having exactly the same amplitude function but time shifted with respect to the top 
y(t) = — t) = where 7 is the input photon bandwidth. This state has zero 

held amplitude, (ain(t)) = 0 , so conventional (second order) interference cannot be used. 
However (aj^(t)ain(t)) 7 ^ 0 so fourth order interference will rehect the quantum coherence 
inherent in the pure state | 1 ^). 

The probability of one and only one click occurring at both detectors Di and D 2 is given 
by the fourth order correlation function 


G(2)(r) 


/o°°/o°°(^lout(^)4,out(^0«2,out(tOai,out(t))dtdF 

/o“(«pout {t)ai,ont{t))dt /“(^)02,out {t))dt 


(4) 


It should be noted that in this expression the time r is not the delay between detection 
events but a temporal separation of the two input photons. In practice, the integration time 
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need not and should not be infinite as it sets the time interval between successive pulses. In 
fact the integration time needs to be of the order Trep ~ max{l/K, I/7}. In what follows we 

h\j 

work in regimes where — > 1, which is compatible with available experimental realizations, 

7 

so we have Xrep ~ I/7. Through equations Q, (|^ and ^ the explicit dependence of 

on g can be seen. By monitoring changes in we can infer changes in g. In the ideal 

case, we would like to detect both photons. 


In practice, either one or two photons could fail to be counted at the detectors due to 
optical losses, imperfect single photon sources or non-unit detection efficiency. The case 
of photon loss before detection (or detector inefficiency) can be modelled as usual [3l] by 
inserting a beam splitter with transmissivity amplitude tj. in the path of the output fields 
The field that reaches an ideal detector is given by the transformation ai^out(^) —t 
tkO-i,out(t) + '>^k'Vi,out(t) where v(t) is a vacuum field mode annihilation operator. Substitution 
into the equation (|^ shows that the only term that contributes to the both the numerator 
and denominator are multiplied by the factor T 1 .T 2 , where = ItfcP is the conditional 

probability for a single photon in the output field to reach an ideal detector. This is because 
this average is normally ordered. Thus (r) is unchanged by loss since we have normalised 
it by the intensity that actually reaches the detector from each mode: in effect G^‘^\t) is 
a conditional probability conditioned on only those detection events that give two counts, 
one at each detector. Single counts and no counts are discarded. These cases should be 
considered as failed, but heralded trials in which we discard and simply run again with 
another two single photons. However, this lowers the sensor’s bandwidth. In section |111|, we 


explicitly include an optical loss factor, e, which is defined as the number of failed trials over 
the total number of trials, to take the effect of these imperfections into account. As photon 
loss is heralded, LOQC error correction techniques might be employed to mitigate the loss 
of signal on such events. 


Photon loss before the device, or failure of a source to produce a photon can likewise 
be modelled by inserting a beam splitter into the path of ain(f). This changes the in¬ 
put state to a mixed state as follows. The input state is a two photon state of the form 
®iin(^)® 2 in(^)| 0 )- W^e then transform ai,in(t) —)■ ai,in(t) = tkai,in{t) + where v{t) is 

a vacuum field mode annihilation operator. Thus the total state after the beam splitter is 
given by aiin(^)h 2 in(f)| 0 ) but the actual input state to the device is given by tracing out 
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over the two vacuum modes. This gives the input state as a mixed state of the form 


Pin — ^1^211)1(11 ® 11)2(11 + 7 "ii? 2 |l)l(l| 0 10)2(01 

+i?iT 2 | 0 )i( 0 | ® 11)2(11 + .Ri.R 2 | 0 )i( 0 | ® 10)2(01. 


(5) 


This also indicates that loss at the input is detected as the conditional input state, condi¬ 
tioned on counting two photons in total at ideal detectors. This is simply the hrst term in 
the above sum. This is the same pure state as for the case of perfect sources. The coefficient 
T 1 T 2 is simply the conditional probability that two input photons in each of the inputs enter 
the device. Input loss or source inefficiency is also heralded in the detectors and those trials 
can be discarded. The fraction discarded in total including input inefficiency and detectors 
inefficiency is simply Ti^inT 2 ,inTi^out 2 ^ 2 ,out, where Tk^in and T^^out are the conditional probabil¬ 
ities that a single photon enters the device in mode-k and that a single photon for output 
mode-k is detected. 

In Fig. the HOM dip for our system is depicted for particular values of K /7 and 
p/ 7 . For r = 0, where input photons are indistinguishable, quantum interference results in 
photon bunching, or photon pairs, and we see the minimum of the coincidence probability 
i.e. the HOM dip. As r increases or decreases the coincidence probability increases. 

We define the responsivity of the sensor to detect the changes in g as 

dG(2)(0) 


Rg{9o,K) = 


dg 


( 6 ) 


Operating at r = 0 is optimal for most combinations of 7 and k and maximizes the 
responsivity. We then optimize the values of K /7 and p /7 so that the derivative of 
with respect to p is maximized. By maximizing the responsivity over our device parameters, 
Po the initial beam splitter coefficient and the cavity damping rate k, we can optimize the 
performance of our sensor. Due to the fact that our sensor is a linear quantum system, we 
can analytically calculate and its derivative for the initial state |'?/( 0 )) = ^ 0 . 2 ,g)-, 

the full expression for G‘^{t) is given in appendix A. Figure]^ can serve as a guide for exper¬ 
imental implementations and device fabrication. Figure |^a) shows as a function of 

p /7 and K/ 7 . Figure [^b) shows the behaviour of the system response for different operating 

points p /7 and K/ 7 . The dashed line on Fig. I^b) demonstrates the operating points at 
dR ^ 

which —= 0 where we can take advantage of maximum sensor response. In addition, 
dg 








FIG. 2. Sensor response to variations in g. (a) Shows the behavior of the estimator, coincidence 
detection probability for indistinguishable input photons versus 5/7 and K/ 7 . (b) Shows 

how responsivity of the sensor varies by operating the sensor at different regimes of g'/y and K/ 7 . 
The white dashed lines show the operating points for which sensor response is maximum and linear 
over the range of small changes in signal. 


at this maximum sensor response, the estimator G^{0) behaves linearly with small signal 
variations as will be described below. 


III. NOISE CHARACTERISTICS 


Another important measnre in characterizing the sensor performance is the Linear Dy¬ 
namic Range (LDR) which is related to the estimation error and the sensor linearity which 
we now explore. The error in estimating 5g is related to the error in estimating (0) in 
a finite nnmber of samples 


5g 


noise 


dG(2)(0) 

dg 


5G'(2)(0)noise, 


(7) 


where 




v/(T(2)(0)(l -G(2)(0)) 

- e) 


( 8 ) 


is the standard deviation of a Bernonlli distribntion with N trials. The loss factor, e, has 
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(a) 



9ly 


(b) 



FIG. 3. Linear dynamic range, (a) Shows how sensor response changes at different operating points. 
If we operate the sensor on a bias where sensor response is maximum, we can take advantage 
of the sensor linear response, up to small variations in g. (b) Shows LDR for bias y = 1.8 shown 
in (a) for different detection frequency bandwidths over 7 , The red star shows the upper LDR 
limit that is the point up to which sensor responds linearly within 1 % variation. 

already been introduced in section [TT| The minimum detectable shift in g from the bias g^ 
should be larger than this error, i.e. ^5'min > <^5'noise, so that we are able to measure it. For a 
large number of samples [N —>■ cx)), ^f/noise is negligible (up to accidental coincidences caused 
by dark counts or stray light). This result is useful for the estimation of a static or quasi 
static parameter. 

We now give an order of magnitude estimate for when the parameter is time varying. 
If ^meas = Rep-^ IS the time between our samples of naive arguments from the Nyquist- 
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Shannon sampling theorem imply that we can not determine freqnency components of g{t) 
greater than / = l/(2Tmeas), which is called the detection freqnency bandwidth. For a one- 
sigma level of conhdence we shonld have N > minjy, k}/( 2/) and the noise eqnivalent 6g, 
given in eqnation ([^, becomes 


^ \/2/GP)(0)(l - GW(0)) 

(^9mm ^ ;— - Yj- -^ 

- e) 


(9) 


Now we can calcnlate the LDR which is dehned as 


LDR = 201og^^, (10) 

^ymin 

where 6g^Q_^ is the point bellow which the sensor response is linear within 1 % variation, i.e. 
Rg{g) = in which = maxRg{gQ, k). 

In Fig. I^a), responsivity is plotted with respect to g for an arbitrary valne of k. We 
bias the initial conpling between the optical resonators (^fo) where the responsivity peaks. 
Therefore, there is a range oi 5g = g — go for which the sensor behaves linearly. LDR is 
shown in Fig. [^b) for some arbitrary detection bandwidth in nnits of 7 . For smaller choices 
of //y, hfy'noise will be decreased, so the sensor can resolve smaller shifts in g. 


IV. HOM SENSOR IMPLEMENTATION 

We now consider specihc physical applications for onr sensor, hrst as a force sensor and 
then as a refractive index sensor, employing conpled L3 PhC cavities [81-IT2] experimental 
data to estimate its responsivity and minimnm resolvable shift in signal for each case. By 
examining the normal mode splitting reported in these references we infer the conpling 
strength g between the PhC resonators is of the order of 10^^ — 10^^ Hz. The evanescent 
conpling strength between the resonators and wavegnides k can be tailored, so that k ~ for 
example. Operating as a force sensor, the measnred signal is the shift in cavity separation 
indnced by an applied force or a strain, while operating as a refractive index sensor, the 
signal to be measnred is a change in refractive index indnced by the presence of a molecnle 
dropped on the air holes between the PhC resonators, for a constant bias cavity separation. 
A shift in either cavity separation, call it x, or refractive index, call it n, modihes the conpling 
strength between the resonators which will be detected by measnring G^(0). Therefore, to 
give an order of magnitnde estimation of the responsivity and minimnm detectable signal 
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in each case we need to investigate the dependence of as a function of x and n. To do 
this we used a ID model analysed by the transfer matrix method [3^ (see Appendix B) 
to investigate the dependence of the cavity normal mode splitting on the change in cavity 
separation or refractive index. 

In the case of identical resonators, Ui = 002 = OJ and ^1 = ^ 2 , the splitting in frequencies 
of the symmetric and asymmetric normal cavity modes is AD = 2g ilES]. Therefore, we 
can write g = ttcAA/A^, where c is the speed of light and A is the cavity mode wavelength. 
Since vrc/A^ is a constant, to hnd the functionality of g with x and n, we need to hnd the 
functionality of AA with those parameters. Numerics show that an exponential function of 
the form g = ae~^^ hts very well on data achieved for normal mode splitting change versus 
different cavity separations (see Appendix B) and an exponential of the form g = 
can describe the changes with respect to refractive index (see Appendix B). Hence, we can 
generally write g{x,n) = We extract the coefficients a,b and d by htting data 

from hgure 2 of citation na for a PhC made of GaAs/AlGaAs (see Appendix G). According 
to their data g is on the order of 10^^ — 10^^ Hz for this range of Xbias that is shown in Fig. 
m We have chosen k of the order of 10^^ Hz. 

A. HOM force sensor 

First we investigate the efficiency of our system operating as a force sensor. In this case 
n = 1 , so by substituting g{x, 1 ) into equation (|^ we can see how the probability of joint 
detections changes for different operating points a;bias (Fig- [^a)). The sensor response 



to changes in x is calculated 


shows that for an 


input photon bandwidth on the order of 7 = IGHz, which is experimentally feasible at the 


moment [371 [38], sensor response to shifts in x is of the order of 10“^(nm)“^. Minimum 
detectable x can be easily related to Sgrain as femin = ^^5'min- Figure ^c) shows this 
noise equivalent x is of the order of 10 “^(nm/A/Hz). Young’s modulus for GaAs is E = 
85.5GPa [3n|- Therefore, for the given lattice with a thickness of t ~ 1pm the stiffness of 
GaAs is A; = y ~ 85.5™. Minimum detectable force is shown in Fig. |^d) and is of the 
order of 10 “^N which compares rather well with the high resolution PhG force sensors [IQIIII] 


exploiting coherent light. However, these schemes use signihcantly larger input power while 
in our results not only is the pulse power (1 ph/pulse) low but also the average power ( 10 “^° 
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FIG. 4. Performance of HOM sensor as a force sensor. This figure is a fabrication guide to building 
a HOM force sensor with maximum performance, (a) Shows how the estimator evolves by changing 
the operating point Xbias for an input photon bandwidth of 7 = IGHz. (b) For the given 7 and k, 
the responsivity is of the order of 10“^(nm)“^. The white dashed lines show the operating points 
where system response to displacement shift is linear, (c) Shows that the minimum detectable 
change in distance is of the order of 10 “^nm/-v/Hz). (d) For PhG made of GaAs/AlGaAs, the 
given value for minimum detectable x corresponds to minimum detectable forces of the order of 
10“^N. Our calculations show that as we reduce kappa, gradually we loose the linear behaviour of 
the sensor (white dashed lines) for smaller Xbias as the best bias points shifts towards larger x or 
smaller g without improving or decreasing the sensor resolution. 
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FIG. 5. Performance of HOM sensor as a refractive index sensor. This figure is a fabrication guide 
to building a HOM refractive index sensor with maximum performance, (a) Shows responsivity of 
the refractive index sensor for different operating points Xbias for an input photon bandwidth of 
7 = IGHz. The white dashed lines show the bias points where sensor response changes linearly for 
very small changes in refractive index. Our theory predicts that responsivity does not depend on 
single photon band width 7 . (b) Predicts that for 7 = IGHz the minimum detectable refractive 
index shift is of the order of 10 “®RIU/\/Hz. 

W), which is defined by the emission rate of the current single photon sources (~ GHz), is 
also low. 

Importantly, fabricating the cavities with a smaller k does not affect the sensor resolution 
but shifts the optimum operating points at which the sensor behaves linearly (white dashed 
lines in Fig. towards larger a;bias- 


B. HOM refractive index sensor 

To operate the system as a refractive index sensor we operate at a fixed x^as, so g{n) = 
_ System response to refractive index shift is Rn{xo,K) = and the 

minimum detectable refractive index shift is calculated as ^rimin = (5fi'min/2(i(7n)|„=i. Figure 
[^a) is a fabrication guide for 7 = IGHz to find the best operating points to achieve max¬ 
imum responsivity together with linear response. Figure |^b) predicts a resolution of the 
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order of 10 ® refractive index unit (RIU) per a/Hz for single photon bandwidth of 7 = IGHz. 
Up to the best of our knowledge the best resolution achieved in schemes Hang is of the 
order of 10“^RIU per \/Hz, however these use more input power. 


V. CONCLUSION 

In conclusion we have described a uniquely quantum protocol for a PhC sensor based 
on two coupled cavities. Our proposal uses single photon states, not coherent states, and 
operates on Hong-Ou-Mandel interference; a fourth order interference effect. The visibility 
for such a HOM proposal is dependent on changes in the coupling between the cavities, which 
is in turn dependent on shifts in the cavity separation distance and/or refractive index of 
the medium in between the two cavities. Very small changes in such parameters result in 
modulation of the cavity coupling rate that can be observed by measuring the resulting 
change in the HOM dip. Our results predict minimum detectable values for refractive 
index and force changes, 10 “® RIU per a/Hz and 10 “^ N, respectively. This estimation 
is based on the parameters obtained form the current experimental implementations of 
coupled L3 PhC cavities in GaAs/AlGaAs [I2] and is not specihcally optimized for our 
sensor. Further development of PhC technology for the sensor presented here could offer 
signihcant improvements in the performance. Our results show that high sensitivity can be 
reached upon achieving high repetition rate single photon sources. 

The advantages of the presented scheme are as follows. This scheme can be implemented 
on chip and fabricated in micro-scale dimensions. Moreover, unlike sensing approaches based 
on transmission spectrum of a L3 cavity coupled to a waveguide, this approach does not 
require spectral resolution that reduces the bandwidth. Additionally, this scheme can be a 
multi-purpose sensor. In this article, we have discussed force (strain) and refractive index 
sensing. With minor modifications, it can be used for other targets such as local temperature, 
pressure and particle detection and analysis. 

A a/2 improvement in estimation accuracy over schemes using coherent light and a \/2 
improvement in bandwidth over schemes using serial single photons can be achieved. The 
improvement of a/2 in estimation accuracy may sound underwhelming. But a different ac¬ 
counting philosophy shows why this factor is important (see Appendix D for further details). 
Suppose we want our estimate of g to have a mean square error of order 10“^ and we ask how 
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many experiments, on average, we must perform to acheive this precision. Then coherent 
light would require 50 x 10® experiments; a serial single photon approach requires 24 x 10® 
experiments; our HOM approach requires 12.5 x 10® experiments. That is we have quartered 
the number of required experiments relative to the coherent light case and halved it relative 
to the serial single photon case. Due to this the reduction in samples the bandwith of our 
sensor relative to both cases is increased. 

The disadvantage of this single-photon-based scheme compared to those using coherent 
light is the difficulty in building reliable single photon sources and detectors. 

To summarise, the key point we are making in the paper is that if one has already 
made a commitment to single photonics in order to gain access to the quantum information 
processing that this provides, single-photon metrology can also be added to the suite of 
tools. HOM interference is the key phenomena that enables scalable quantum information 
processing in single photonics with linear optics. 


VI. APPENDIX A 


The solutions to the quantum Langevin equations (|^ are given by 
ai{t) = ^/K A{t) [ ciT(c'(T)ai,in(t') + D{t')a 2 ,in{t') 


+B{t) / dt'(D{t')ai;in{t') + C{t')a2;m{t') 


aiit) = a/k Bit) / dt'+ D{t')a2^init') 


+A{t) J dtyD{t')ai^iriit') + C{t')a2^mit')j , (11) 

where A{t) = cos{gt), B{t) = sui^gi), C{t) = cos{gt) and D{t) = 

^gKt/2 sin(^^)_ By using the above solutions for the cavity mode in the input-output relation 
(|^, we can analytically calculate the joint detection probability as 


--t(k+7) / 3 1 1 

/o\ € ^ I —r(K+7) — —T(3K-\-'y) —r(«:+37) / , \ i 

G(2)(r) =- j - I He2 ’ + Ce 2 ^ +De2 ’ + | , (12) 


where 


A = (4/ + + (f - + ^g\i^ + «:^) 1 , 
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B = { Ag ^ + (7 — ^ 256 ( 7 ® + ^^(7 + kY + §(7^(7^ — 2 k ^)( 16 ( 7 '^ + ^^(7 + kY ) 

+ 16/(7^ + 27^^^ + 207^3 + 22 k ^)^ , 

C = -32/k2(V + 7^ - 

D = -32c/V«^^^^ 

E = — 645(^7k^(4(7^ + 7 ^ — k^)( 4(7^ + /s:^)F, 

and 

F = /?(—12(7^ — 7^ + cos(77r) + 2 g{Ag'^ + 7^ — 3 k^) sin((7r). 

VII. APPENDIX B 

To find the fnnctionality of conpling strength g with separation distance between the 
cavities and refractive index of the media in between the two cavities, we can nse the 
analogy of the conpled cavities with a qnantnm donble-well problem with a potential barrier 
in between, where we need to hnd the tnnnelling rate g. Solving this problem shows the 
fnnctionality of g with the width of the barrier x and the height of the barrier V scales as 
g cc exp{——dV}. According to citations [lUllS], which introdnce the optical eqnivalence 
of the Scrodinger eqnation, V oc —n^. Therefore, we expect the conpling strength between 
the cavities to scale as g (x exp{—6a: + To be more precise a one dimensional optical 

modelling simulation has been performed to hnd the functionality of coupling strength g with 
X and n, which suggest a very good ht can be achieved with the above given functionality. 
The simulations are done by the transfer matrix method described in [35]. The results are 
shown in Fig. 
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FIG. 6. Approximating the optical coupling between two GaAs cavities (L = 450 nm) placed in the 
middle of a distributed Bragg reflector (DBR) stack comprising GaAs (d=225 nm) and air (d=112 
nm) pairs versus cavities separation (a,c) and air holes refractive index (b,d). (a) Normalised 
transmission of the stack showing the normal modes of the two cavities for different x. The dashed 
line corresponds to the unperturbed single cavity mode confined by the DBR stack, (b) Normal 
modes shown for x=449 nm (two air and one GaAs layers) and varying the refractive index of the 
air layers between the cavities. By increasing the air hole refractive index the mode separation 
increases which corresponds to stronger coupling between the cavities due to the decrease in the 
refractive index offset of the DBR layers, (c) coupling frequency calculated from (a) versus cavity 
separation. The dashed line corresponds to an exponential decay fitting, (d) coupling frequency 
as a function of air hole refractive index calculated from figure (b). The dashed curve corresponds 
to an exponential fitting of . 
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VIII. APPENDIX C 


We choose experimental data for a PhC lattice made of GaAs/AlGaAs given in [T2] as 
an example to find the functionality of coupling strength g in terms of cavities separation 
X and the refractive index of the dielectric material n. Figure shows the numbers for the 
best found htted function. 



FIG. 7. Functionality estimation of g versus x and n. Fit on experimental data (red points) given 
in figure 2 of citation |12j to find the coefficient a, b and d in functionality of g versus x and n that 
we found of the form g{x,n) = vrcAA/A^ = where A = lOOOnm. 


IX. APPENDIX D 

Here we back up the claims made in the conclusion, further details will be available in 
a future publication. Here we assume a mean flux greater than two photons may damage 
a sample. We model our coupled cavity sensor as an effective beam splitter between the 
input modes such that the beam splitter reflectivity g is a. function of the parameter to be 
estimated, x, i.e g{x). In the case of very small changes in x, this functional relationship 
can be linearized about an operating point xo so that changes in 6 are linearly proportional 
to changes in g. 

Now we ask what is the ultimate limit imposed by quantum theory on the precision of our 
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estimate gest of g. This, of course depends on the state used to probe the beam-splitter. If 
the mean squared error (MSE) i.e. E[(5f —quantifies the performance of our estimation 
scheme then the quantum Cramer-Rao bound provides a method to answer this question. 
The quantum Cramer-Rao bound states that asymptotically the precision of an unbiased 
estimation scheme is bounded below by 1/■>/ NI{g) where N is the number of experimental 
trials and I{g) is the quantum Fisher information with respect to the input state. 

If the input state is a coherent state then the quantum Fisher information is generally 
Ia{g) = where |ap is the mean photon number at the input. So for the cases of 
interest Ia{g) = 1 and Ia{g) = 4. If the input state is a Fock state then the quantum Fisher 
information is [M]: I pig) =4 (input photon number = 1), Ip{g) = 16 (input photon number 
= 2). The numbers quoted in the main text can be arrived at by setting MSE = 10“^ and 
solving for N. 

Loss is easily included using the standard beam splitter model discussed in section 2. 
Coherent states are not entangled on beam splitters so tracing out the vacuum modes leaves 
the signal mode in a coherent state with amplitude reduced by a —)■ ta, where = T is the 
probability for a photon not to be lost. The Fisher information |16] is then simply rescaled 
to reflect the loss of amplitude. Single photon product states at the input to a beam splitter 
do become entangled at the output. However, we can trace out the vacuum modes to give 
a mixed input state for the case of n = 2 in a HOM experiment of the form 

Pin = ^ 1 ^ 211 ) 1(11 ® 11 ) 2(11 + ^i-R2|1)i(1| ® 10 ) 2(01 (13) 

+ -RiT 2 | 0 )i( 0 | ® 11)2(11 + -Ri-R 2 | 0 )i( 0 | ® 10)2(01. 

The Fisher information for a mixed state is more difficult to calculate as it involves the 
logarithmic derivative. However, because the mixed state above is so simple it is easy to 
see that, conditioned on detecting two photons for correct heralded operation, the Fisher 
information is rescaled by TfT^. 
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